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ABSTRACT 

We present a detailed analysis of the baryonic and dark matter distribution in the lensing cluster 
Abell 611 (z=0.288), with the goal of determining the dark matter profile over an unprecedented 
range of cluster-centric distance. By combining three complementary probes of the mass distribution, 
weak lensing from multi-color Subaru imaging, strong lensing constraints based on the identification 
of multiply- imaged sources in Hubble Space Telescope images, and resolved stellar velocity dispersion 
measures for the brightest cluster galaxy (BCG) secured using the Keck telescope, we extend the 
methodology for separating the dark and baryonic mass components introduced by Sand et al. (2008). 
Our resulting dark matter profile samples the cluster from ~3 kpc to 3.25 Mpc, thereby providing an 
excellent basis for comparisons with recent numerical models. We demonstrate that only by combining 
our three observational techniques can degeneracies in constraining the form of the dark matter profile 
be broken on scales crucial for detailed comparisons with numerical simulations. Our analysis reveals 
that a simple Navarro, Frenk, and White (NEW) profile is an unacceptable fit to our data. We confirm 
earlier claims based on less extensive analyses of other clusters that the inner profile of the dark matter 
profile deviates significantly from the NEW form and find a inner logarithmic slope /3 flatter than 0.3 
(68%; where pdm oc r^^ at small radii). In order to reconcile our data with cluster formation in a 
ACDM cosmology, we speculate that it may be necessary to revise our understanding of the nature of 
baryon-dark matter interactions in cluster cores. Comprehensive weak and strong lensing data, when 
coupled with kinematic information on the brightest cluster galaxy, can readily be applied to a larger 
sample of clusters to test the universality of these results. 

Subject headings: dark matter — galaxies: clusters: individual (Abell 611) — galaxies: elliptical and 
lenticular, cD — galaxies: formation — gravitational lensing 



1. INTRODUCTION 

The cold dark matter (CDM) model has been remark- 
ably successful in explaining the observed large-scale 
structure of the universe (e.g. JSpringcl ct al. 2006). Cos- 
mological A^-body simulations have been a crucial tool in 
assessing CDM, making precise predictions of the growth 
of cosmic structure and the distribution of dark matter 
(DM) over a wide range of scales. The density profile 
of DM has particularly been the focus of intense theo- 
retical and observational efforts, as it offers a tractable 
route toward testing the nature of dark matter and its 
interaction with baryons. Key issues include the form of 
the density profile, its universality over a wide range of 
halo masses, the degree of variance among halos, and the 
role of baryons in shaping the DM distributi on. 

Ba sed on A^-body simulations, Navarro ct al.l ()1996l . 
IT997I hereafter NEW) found that CDM halos are self- 
similar, differing only by simple rescalings of size and 
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density, over 4 decades in mass. Within a scale radius 
Ts, this "universal" NEW density profile asymptotically 
approaches p oc , with /3 = 1; on larger scales, it 
approaches p oc . Moreover, NEW found to be 
correlated with mass, implying that the diversity of DM 
halos is captured by a single parameter. Subsequent 
work with higher num erical resolution con firmed t he ba- 
sic fi i idings of NEW (iMoore et a l.l 119981 : Ghi gna et all 
120001: iDiemand et al.l l2004. 2005). but disagreed on the 
inner logarithmic slope (3. Nonetheless, a steep {(3 > 1) 
cusp was favored by all authors, with claims ranging up 
to /3 « 1.5. Subsequently, a generahzed NEW (gNEW) 
profile was introduced: 

P^''^ ^ (r/r,)/5(l + r/r,)3-/'' 

This reduces to the NEW form when /3 = 1. More re- 
cent studies by iNavarro et~an ()2004l . I2008f l have hinted 
that, rather than converging to a power law, the den- 
sity slope may become progressively shallower at very 
small radii < 10~^rs. These authors also emphasized 
slight differences in form among halos, demonstrating 
that self-similarity may not strictly be the case; how- 
ever, it must be emphasized that these deviations from 
NEW, although significantly determined, are still small. 

Despite its triumphs on cosmological scales, CDM pre- 
dictions on galaxy scales have often been difficult to rec- 
oncile with observations. A particular challenge is the 
high central density expected in DM halos. Studies of 
dwarf and low surface brightness galaxies, which are ex- 
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pected to be DM dominated throughou t, reveal kinemat 
ics incompatible with CDM hales (e . g..lCote et al.|[2000l: 
Marchesini et al.ll2002t LS imon e t all l2003t iGentile et al ' 



2005; Kuzio d e Naray et al. 200^ but see IS waters et al 
2003) . These data are often more consistent with a shal- 
low cusp (/3 < 1) or cored (/3 = 0) density profile than 
with an NFW-like cusp. Rotation curves of spiral galax- 
ies have also indicated lowe r central den sities than are 
expected in CDM halos (e.g..lNavarro fc S tcinmctz 2000; 
Mo & Mao' '2003; iSaluccil I2001^ 'Binncv & Evana i2001l : 
Alam et al. 2002; McGaugh et al. 2007). It is generally 
thought that accounting for baryons will make simulated 
halos even more concentrated; we discuss the current lit- 
erature in fj6l including contrary indications. 

Thorough observational testing of the numerical sim- 
ulations requires checking their predictions on all scales. 
Galaxy clusters provide an excellent laboratory for 
these tests because several independent observational 
probes are available, spanning a large dynamic range 
in density and cluster-centric radius. The mass pro- 
file of clusters has been a hotly debated topic in re- 
cent years. The inner slope (3 and the concentration 
parameter c = rmr/fs have been particularly contro- 
versial, with different groups obtaining results ranging 
from perfect agreement to highly significant disagree- 
ment between observations and CD M predictions, of- 
ten for the very same clusters (e.g.. Smith et al. '20051 
'^avazzi 2005; Zappacosta ct al. 2006; Schmidt & Allep 
2007; Umetsu & Broadhurstl 120081: IBradac et all 120081 
Limousin et al..i2008l ). 



In addition to scatter in the mass profile from cluster to 
cluster and possible selection biases, other factors clearly 
contribute to this broad and often conflicting range of re- 
sults. First, different groups adopt different definitions 
for the inner slope /3; some use it to refer to the total 
mass density profile, whereas others use it to refer only 
to that of the dark matter as we do in this paper. Second, 
the stellar mass of the central galaxy is often neglected 
on the grounds that it makes a negligible contribution to 
the total mass. However this is not the case in the inner- 
most regions where stars typically dominate the density 
profile. Third, studies that rely on an individual observa- 
tional technique cover a limited dynamic range in radius, 
thus effectively relying on an NFW (or other) model to 
extrapolate the results well beyond the region probed by 
the data. 

X-ray temperat ure measurements typically reach out 
to ~500 kpc fe.g.. lBradac et al.ll2008f) . and are often lim- 
ited to radii larger than ~50 kpc due to instrumental res- 
olutio n or residual substructure (e.g., iSchmidt &: AllenI 
l2007f ). Strong lensing is typically sensitive to the pro- 
jected mass distribution inside '^100-200 kpc, with lim- 
its on at '--^10-20 kpc available in only the most fa - 
vorable cases (e.g., ,Gavazzl,2005t ILimousin et al.ll2008n . 
Weak lensing requires averaging noisy signal from many 
background galaxies and therefore does not have the 
resolution to constrain profiles inside ~100 kpc. Stel- 
lar kinematics of the central galaxy can cover the 
^1-200 kpc region with l ong exposures on large tele- 
scopes ([Kelson et al. I [20021 ). Satellite kinematics is typ- 
ically limited to rad i i larger than ~ 100-200 kpc (e.g., 
iDiaferio et all 120051 : iBiviano fc Saluccil I20060 . It is 
thus clear that only by combining multiple diag nostics 
(e.g, lMiralda-Escudc.1995. : .Sand et al...2002. : .Kneib et all 



[200l [Bradac et al.l [20051 : IMahdavi et~all[200l can one 
hope to achieve a precise and accurate determination of 
the inner slope, and thus assess the validity of the NFW 
profile. ^ ^ 

W this goal in mind, ISand et al.l (|2002l . 120041 ) stud- 
ied a sample of nearly round, apparently "relaxed" 
clusters by combining strong gravitational lensing con- 
straints with the velocity dispersion profile of the cen- 
tral brightest cluster galaxy (BCG). The combination 
of the two techniques allowed them to disentangle the 
contribution of the stellar mass of the BCG from that 
of the dark matter halo and to obtain constraints over 
almost two decades in radius. They found that the 
steep cusps produced in A^-body simulations were in- 
consistent with the data, and instead favored a shal- 
lower slope with (/?) — 0.52 ± 0.05 ± 0.2 (statistical and 
systematic errors). Following claims by some workers 
that this discrepancy could be an artefact of simplifying 
assu mptions in the analysis, such as negligible elliptic- 
itv (iBartelmann fc Meneghettil 120041 : iMeneghetti et all 
120071) . ISand' et al.l (|2008l ) presented a more sophisticated, 
fully two-dimensional analysis of two clusters. They 
reached similar conclusions (e.g., (3 = 0.451q25 for Abell 
383), but again noted degeneracies (e.g., between and 
(3) that could only be broken by including mass tracers 
at larger radii. This paper is th e next step in a continu- 
ing effort to build on the work of Sand et all ([20021 . 12004 
[200 8.V We are currently collecting data for a sample of 
10 clusters, and present here a p ilot study of Abell 611 
{z = 0.288. iCrawford et al.|[T995f) . This cluster exhibits 
a regular, nearly round morphology in X-ray and opti- 
cal imaging. In addition to the strong lensing constaints 
(30-90 kpc) and stellar kinematics (< 20 kpc), we have 
obtained wide-field Subaru imaging to measure mass 
on 0.15-3 Mpc scales via weak len sing, thus b reaking 
the remaining degeneracies noted in iSand et all HOOS). 
The goals of this paper are to determine whether an 
NFW profile provides an acceptable fit to all constraints, 
whether the gNFW profile substantially improves the fit, 
and if so, the range of inner slopes (3 that is permitted. 

A full comparison of observational data with numerical 
simulations is possible only by probing mass on all scales. 
It is much easier for theory to fit data acceptably in a 
single regime (e.g., on weak or strong lensing scales alone) 
than over the full extent of the cluster. We thus consider 
the wide dynamic range of our data (3 decades in radius) 
to be a prime advantage of our analysis. Additionally, a 
proper comparison with theory requires separating DM 
from baryons. Although baryons are a small fraction of 
the total cluster mass, they can dominate on kpc scales. 
We therefore carefully model starlight in the BCG to 
account properly for stellar mass. 

The paper is organized as follows. The next three sec- 
tions (§[2|to[l|) describe the three main datasets used in 
the analysis. Each dataset is the foundation of one of 
three mass probes that are combined in this paper. Dis- 
cussing each one separately allows us to present the mea- 
surements in detail and the strengths and weaknesses of 
the probe: 321 presents a weak-lensing analysis based on 
wide-field Subaru imaging; iJ3| presents a strong lensing 
analysis based on Hubble Space Telescope data, which 
is also used to constrain the distribution of starlight on 
galactic scales; |j4| presents a dynamical analysis of the 
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stellar kinematics in the BCG obtained from Keck spec- 
troscopy. Our results on the dark matter density pro- 
file of the cluster based on the combination of the three 
probes is presented in fJSJ illustrating how this power- 
ful combination breaks the degeneracies inherent to each 
method alone. Limitations of our analysis and residual 
systematic effects are also considered. SJSl discusses our 
results in the context of the literature, and SJT] summa- 
rizes. 

Throughout the paper we adopt the cosmological 
parameters (/i, r2,„, il^) — (0.7,0.3,0.7). Magnitudes 
are given in the AB system unless otherwise stated. 
Marginalized posterior distributions are characterized by 
the mode, and error bars refer to the 68% confidence 
level. 

2. SUBARU IMAGING AND WEAK LENSING 
ANALYSIS 

This section presents an analysis of wide-field, multi- 
color Subaru imaging data whose purpose is to measure, 
via the gravitational shear, the radial mass distribution 
from 150 kpc to 3.25 Mpc. The shear signal is limited 
outside this interval by the surface density of suitable 
background galaxies on small scales, and by confusion 
with large-scale structure on large scales. We use photo- 
metric redshifts to identify foreground and cluster mem- 
ber galaxies, which are not lensed by the cluster and so 
dilute the shear signal if not carefully excluded. Shapes 
of galaxies are measured and corrected for atmospheric 
and instrumental distortions using the method originally 
outlined by Kaiser ct al. (199^, hereafter KSB). We have 
calibrated this procedure using simulated data. The two- 
dimensional mass map reveals no significant secondary 
mass concentrations, justifying our use of simple para- 
metric models. We show that our radial mass profile is 
in excellent agreement with independent X-ray observa- 
tions. 

2.1. Observations and Data Reduction 

We observed Abell 611 on 2007 November 12-13 us- 
ing SuprimeCam at the Subaru Telescope under excel- 
lent conditions. BjVjRcIc imaging was obtained, with 
1200 s integration in Vj and 2400 s in each of BjRcIc- 
We conduct shape measurement in the the Rc imaging; 
the remaining filters are used for photometric redshifts. 
The median stellar FWHM in the Rc imaging is 0'.'7, 
and the la surface brightness limit is 27.2 mag arcsec"^. 

The SuprimeGam imaging was reduced using the Im- 
CAT^ software. One corner GGD (DET-ID 0) was ex- 
cluded from our analysis due to poor charge transfer ef- 
ficiency. Bias was measured from the median of 6 frames 
and subtracted. Saturated pixels, chip defects, ghosts, 
and satellite trails were masked, and cosmic rays were 
rejected using gethotpix. Since scattered light is sig- 
nificant and non-uniform across the camera, flat-fielding 
must be conducted on a chip-by-chip basis. In each fil- 
ter, a flat was constructed from the night sky imaging 
by placing large elliptical masks around objects detected 
by hf indpeaks, normalizing each frame by its median, 
and then median-combining all frames with a 2a clip. 
In order to have enough frames for this procedure to be 

^ |http: //www, if a.hawaii . edu/~kalser/imcat/| 



effective, we combined data from both nights. The chip- 
to-chip normalization was provided by the SuprimeGam 
team (H. Furusawa and Y. Komiyama, private commu- 
nication) . 

An initial sky subtraction is made by subtracting the 
mode of each chip. However, large-scale gradients remain 
from residual scattered light and Galactic cirrus. We 
therefore adopt a more sophisticated backgro und sub- 
tractio n sc heme adapted from that described bv lDonovanI 
(|2007| ) and lGapak et all (|2007D . Roughly, each chip is di- 
vided into a grid, and the background is measured in each 
cell after excluding objects. The cells are interpolated via 
Delaunay tesselation to form a continuous background 
image, which is then subtracted. This method is suc- 
cessful at removing background gradients, stellar halos, 
and charge buildup. 

Photometry was calibrated using two Landolt 
fields. Herve Aussel kindly provided AB magni- 
tudes in the SuprimeGam filters for our fields, us- 
ing fits to PH OENIX stellar atmo sphere models 
(iHauschildt et al.l IT 997) with U BVRI (jLandoltl [1992). 
ugriz (Adelman-McGarthv ct al."2008', SDSS DR6), and 
JHK (Skrutskic ct al. 2006, 2MASS) photometry. We 
used 5 stars with well-fit SEDs to measure the photo- 
metric zero-point, and typical Mauna Kea values for the 
airmass correction, which was always < 0.03 mag. 

The astrometric solution was obtained by fitting the 
positions of unsaturated stars to a low-order polynomial 
with iterative rejection. The median dispersion among 
stellar positions in different exposures was 0'.'007. Ab- 
solute astrometry (tangent point, rotatio n, scale) was 
ancho red to stars in the USNO-B catalog (|Monet et al.l 
[200l . Ima ges were placed on a stereographic projection 
and combined using a clipped weighted mean. 

Noise maps were created by measuring the variance in 
each chip (before applying distortion) and propagating 
through the coaddition. Diffraction spikes, halos and 
rings of very bright stars {V < 10), and perimeter regions 
were masked in the stacked image. 

2.2. Source Catalog and Photometric Redshifts 

Object detection was performed using SExtractor 
(jBertin fc ArnoutsI 119961 ) with a threshold of 0. oa over 
> 15 contiguous pixels. Extinction was corrected i n each 
filter using the reddening map of ISchlegel et"all ([1991 ) 
and the Galactic dust law. Golors were measured in 
3" apertures, with aperture corrections determined us- 
ing bright starS; 

The BPZ (,Benitezll2OOO0 code was used to obtain photo- 
metric redshifts using the i^^-band MAG_AUTO mag- 
nitude and the color measurements described above. The 
CWWSB templates and a prior on the luminosity func- 
tion from the HDF were used. For our weak lensing 
galaxy sample, cuts are placed on Zf,, the marginalized 
redshift, and on P{z > 0.4), the probability that z > 0.4. 
We require 0.4 < Zb < 1.3, with the upper limit being 
roughly the redshift at which the 4000 A feature redshifts 
outside the Ic filter and beyond which there is significant 
degeneracy with low-redshift galaxies. We also require 
P{z > 0.4) > 0.75; this cut was determined from Monte 
Carlo simulations by requiring 90% confidence in < 10% 
contamination by unlensed galaxies in the inner cluster 
(< 250 kpc projected radius). 

We have checked the redshift distribution we obtain in 
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Fig. 1. — Redshift probability density for R < 25.8 (~ Scr) 
galaxies in Abell 611, as well as 4 comparision fields also ob- 
tained with Subaru, but with wider wavelength coverage: A370 
(UBVRIZ), HDF-N (UBVRIZ, HK, NB816), Lockman Hole 
Northwest {BVRIZ.HK), SSA22 (UBVRIZ JHK). Note spikes 
in the A611 {z = 0.3) and A370 \z = 0.4) fields at the cluster 
redshifts. Our redshift distribution is generally bracketed by the 
comparison fields and appears consistent within the limits of cos- 
mic variance. Galaxies used in our weak lensing analysis lie in the 
redshift interval indicated by the arrow; however, the additional 
cut P{z > 0.4) > 0.75 eliminates most galaxies with z < 0.5. 

Abell 611 against those obtained in other Subaru fields 
with wider wavelength coverage, as shown in Figure [TJ 
The redshift distribution for Abell 611 is bracketed by 
those of the four comparison fields, with no obvious sys- 
tematic errors within the limits imposed by cosmic vari- 
ance. For each field, we have also calculated the effective 
redshift Ze, at which the lensing efficiency Dis/ Dg equals 
the sample mean, restricting to 0.5 < z < 1.3. (The 
lower limit is to avoid the cluster in Abell 370; the up- 
per limit is that which we adopted in Abell 611.) For 
Abell 611 Ze = 0.79, in good agreement with the other 
fields: Ze = 0.78 (A370), z^ = 0.82 (HDF), z^ = 0.82 
(LH-NW), Ze = 0.89 (SSA22). We see no evidence that 
our BVRI coverage is insufficient to provide redshifts of 
sufficient accuracy in the restricted range z < 1.3. Since 
our data and the reference fields have been analyzed sim- 
ilarly, in particular using the same bpz code, they are 
potentially subject to common systematic errors. Other 
studies have found that BPZ performs well in compari- 
son t o large spectroscopic samples (e.g.. iMobasher et al.l 

UMl). 

We conclude Zf. is unlikely to be in error by more than 
0.03, corresponding to an error of 2% in the shear normal- 
ization. The propagation of this uncertainty is discussed 
in §3231 and [Ml 



2.3. Galaxy Selection and Shear Measurement 

Shapes of galaxies must be corrected for both distor- 
tions introduced by the point spread function (PSF), 
which is generally anisotropic and varies across the fo- 
cal plane, as well as the isotropic smearing introduced 
both by the atmosphere and by the smoothing kernel 
used in our shape measurements. The determination of 
the PSF from stellar images and its application to cor- 
recting galaxy shapes is described in Appendix A. The 
interested reader can also find details of our measure- 
ment of the reduced shear ga in Appendix B. Here, we 
describe only the weak lensing sample selection. 



Proper measurement of the weak gravitational shear 
requires careful selection of the lensed galaxy sample, 
avoiding contamination by cluster members and fore- 
ground objects which could dilute the shear signal. We 
adopt the following cuts to remove spurious detections 
and select galaxies that are resolved, significantly de- 
tected, uncontaminated by neighbors, and located well 
behind the cluster: 

• l.l(r^) < r/i < 6 pixels, where r/j is the half-light 
radius and (r,*) is the median stellar half-light ra- 
dius. The lower limit selects resolved objects; the 
upper limit removes foreground objects and spuri- 
ous detections. 

• Positive flux detected with S/N = 
FLUX_AUTO/FLUXERR_AUTO > 5, which tests on 
STEP simulations ( ij2.4p showed to be a reason- 
able threshold. 

• Rc > 21 mag, approximately the saturation limit. 

• PSF-corrected ellipticity |e'| < 1, since in the weak- 
lensing regime, higher ellipticities often indicate 
blends. 

• No masked pixels, or pixels assigned by SEXTRAC- 
TOR to nearby neighbors, within 2rg (defined in the 
Appendix) . 



tr Psh > and tr Psm 
pendix. 



> 0, as defined in the Ap- 



• d < 0.2 pixels, where d is the distance between 
the centroids measured with and without smooth- 
ing. For symmetric objects, d = in the absence 
of noise. We find that large separations usually 
indicate blended or asymmetric objects. 

• 0.4 < z < 1.3 and P{z > 0.4) > 0.75, as discussed 
in go 

This selection yields a galaxy density of 14.5 arcmin"^, 
comp arable to other stud ies with the Subaru telescope 
(e.g.. lUmetsu et al.|[2009t ). 

2.4. Shear Calibration 

The KSB metho d is known to und erestimate shear 
by ~10-15% (e.g., lErben et all [200l . We therefore 
calibrated our met hod using the STEP2 simulation 
(|Massev et al.|[2007l ). in which shapelet models of galax- 
ies are subjected to an anisotropic PSF and known, con- 
stant shear and placed in an image designed to mimic 
Subaru data. PSF models A and C have sizes ~0'.'6 and 
~0'.'8, respectively, which bracket the 0'.'7 seeing in our 
data. We measured the shear in 128 simulated fields (64 
per PSF model) using the same method applied to the 
real data, with a few exceptions. First, noise was esti- 
mated by SExTRACTOR from the variance in the coadded 
image and corrected upward to account for noise correla- 
tion. Also, orders of the interpolating polynomials were 
reduced owning to the smaller 7' field containing fewer 
galaxies. 

For each field, the weighted mean shear and its uncer- 
tainty were estimated and are shown as a function of the 
true shear for PSF A in Figure [2l (STEP2 simulations 
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Fig. 2. — Response of our shear measurement technique to simu- 
lated shear fields from STEP2 PSF A (0^6). Crosses and diamonds 
denote the gi and g2 components, which are best fit by the solid 
and dashed lines, respectively. A typical error bar is shown. 



TABLE 1 
Shear, Calibration With STEP2 
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Fig. 3. — Smoothed galaxy light distribution with shear field and 
weak lensing mass contours overlaid. Contours are linear in mass, 
with the lowest corresponding to ^2(7. Light regions represent 
those around masked bright stars. Within the uncertainty from 
binning, the mass and light centroids are coincident. The width 
of the ICF used in the mass reconstruction is shown; the light is 
smoothed to the same scale. 



PSF 


Component 


c 


m 


A 


91 


(-1.1 ±0.7) X 10-3 


0.86 ±0.03 


A 


92 


(0.7 ±0.7) X 10-3 


0.79 ± 0.03 


C 


91 


(2.3 ± 0.7) X 10-3 


0.80 ±0.03 


C 


92 


(0.5 ± 0.7) X 10-3 


0.80 ±0.03 




Average 




0.81 ±0.03 



Note. — Errors are obtained from least-squares linear 
fitting. The rms residual in each case is 0.004. 



have zero convergence, so there is no distinction between 
shear and reduced shear.) A hnear relation 

^?a,mcas — ^^7a, input C (2) 

was fit to both shear components of each PSF model. We 
note no signs of non-linear response in this shear regime. 
The fit parameters m and c are given in Table[T] Additive 
error reflected in c results from error in the PSF determi- 
nation and is < 10"'^; we conclude that our correction for 
PSF anisotropy is accurate. Multiplicative calibration er- 
ror is reflected in m. The mean and standard deviations 
of these measurements is to = 0.81 ±0.03, and we cor- 
rect all shear measurements and associated uncertainties 
in this work by this factor. This is within the range o f 
m reported for STEP2 participants (|Massev et al.ll2007D . 
although rather at the low end. 

2.5. The Mass Distribution from 150 kpc to 3.25 Mpc 

We use the Lensent2 code (jMarshall et all 120021 ) to 
produce a mass map and verify that the dark matter dis- 
tribution in Abell 611 is apparently relaxed, symmetric, 
and centered on the BCG. A Gaussian intrinsic corre- 
lation function (ICF) with a FWHM of 170" was used, 
since this maximises the detection significance; we also 
note that larger ICFs do not significantly increase the 
evidence. Figure [3] shows the smoothed light distribu- 
tion and mass contours. The offset between the light 
and mass centers is within the bin size used in the mass 
reconstruction, and we thus consider it insignificant. In 
addition to the cluster core, a few other mass concentra- 
tions are marginally detected. The most significant of 
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Fig. 4. — Top: Radial profile of tangential reduced shear and 
the best-fitting NEW model (Table 2). Bottom: S-mode signal, 
consistent with zero (x^/dof = 6646/6545 = 1.02). 

these is approximately 1.5 Mpc north. The excess mass 
in this clump is < lO^'' M0 , or < 15% of the cluster virial 
mass, and the detection is < 3cr. Our treatment of this 
clump is discussed further below. 

Having established that the weak lensing signal can be 
described by a single concentration of mass, we now pro- 
ceed to an analysis based on simply parametrized models. 
For this purpose we consider only the signal arising from 
background galaxies whose projected distance from the 
cluster center is between 150 kpc and 3.25 Mpc. The 
outer limit is motivated by the onset of contamination 
from large-scale structure (|Hoekstrall2003[ ). We adopt the 
inner limit to avoid the strong-shear regime, unprobed by 
STEP2, where non-linear calibration errors may become 
significant; in any case, the number of selected galaxies 
within 150 kpc is very small. The individual photometric 
redshifts of sources are used in our shear calculations. 

Throughout this pape r, we use a modified version 
of the Lenstool code^ (jJullo et al.l l2007t iKneib et all 

^ |http : //www ■ oamp . f r/cosmology/lenstool/ 1 
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TABLE 2 
NFW Fit to Weak Lensing Data 



Parameter 


Units 


Prior 


Posterior 


e 




[0,0.3] 


< 0.16 


PA 


deg 


[0,180] 


28 ± 26 


Ts 


kpc 


[50, 800] 


301+189 


fO.DM 


km 


[1000, 2200] 


1525 ± 110 


M200 


10" M0 






C200 








m 




(At, cr) = (0.81,0.05) 


0.80 ±0.06 


Minimum 


^WL (13090 constraints) 


13083.2 



Note. — Brackets specify the range of uniform priors. 
For the shear cahbration m (Eqn. [2]l, a Gaussian prior was 
used. fTo,DM is defined in Eqn. |4] e is the pseudoellipticity 
of lGolse fc Kneibl pOOl ). 



119931) to fit data to parametrized models. Samples are 
drawn from 10 parallel chains using the Markov Chain 
Monte Carlo (MCMC) technique with simulated anneal- 
ing. The rate parameter, which controls the convergence 
speed, was set to 0.05 in order to encourage full explo- 
ration of the parameter space, and to obtain sufficiently 
precise measures of the evidence with reasonable compu- 
tational effort. In all cases, we repeat the MCMC anal- 
ysis 4 times to estimate the error on the Bayesian evi- 
dence. Parameters were estimated by combining chains, 
using at least 6000 samples. 

For comparison with the literature and with theory, 
we fit the shear data to a pseudoelliptical NFW profile. 
Since the mass of cluster galaxies is negligible on the 
large scales at which we measure shear, the mass model 
includes only the cluster-scale DM halo. The shear cal- 
ibration m is included as a free parameter, with a prior 
based on the STEP2 results of H2.41 thereby propagating 
the calibration uncertainty. Recognizing the potential 
for small systematic redshift errors that mimic shear cal- 
ibration error (i i2.2p . as well as the systematic uncertain- 
ties inherent in applying simulation-based calibrations to 
real data, we have slightly inflated the error in m to 0.05. 
For each model suggested by the MCMC sampler, the re- 
duced shear polar gi, model is computed at the position of 
each galaxy i, and 



E 



Si , mcas S' 



-i, mo del 



(Ti/m 



(3) 



is calculated, where ai is the shear uncertainty described 
in Appendix IB] 

The Lenstool code parametrizes the NFW profile 
with the canonical scale radius and (To, dm, a char- 
acteristic velocity dispersion defined by 



(4) 



The results of the fit are presented in Table [2] and ex- 
pressed both in terms of Lenstool parameters and in 
the standard parametrization of the NFW profile using 
a-iid C2oo-^*' The measured radial shear profile is 
compared with the results of the fit in Figure H) The 
i3-mode signal - consistent with zero - is also shown, 

When elliptical models are considered, we use a circularized 
measure for the radii r200 and r^, and hence C200 = ^200 /r a- 
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Fig. 5. — Radial surface brightness profile of the Abell 611 
BCG with PSF-convolved R^/^ and dPIE fits. The mean sur- 
face brightness in radial bins is denoted with crosses, and the la 
dispersion in each data bin is shaded. The R^^'^ and dPIE fits 
are used in modeling the strong lensing and velocity dispersion 
constraints, respectively. Here the elliptical radius is defined by 
r2 = x2/(l + e)2 + j/2/(l - e)2, where e = (a^ _ 62)/(a2 yiy 



indicating that systematic errors are smaller than statis- 
tical errors. 

We have experimented with a two-component mass 
model incorporating the marginally-detected clump pre- 
viously noted in our two-dimensional mass reconstruc- 
tion. However, Xwl '^ot improved, and the Bayesian 
evidence is much reduced, strongly arguing against the 
need for modeling the clump. We therefore do not con- 
sider any additional large-scale mass components for the 
remainder of this paper. We also note that assuming ax- 
ial symmetry, rather than the two-dimensional elliptical 
symmetry we adopt here for consistency, has a negligible 
effect on the inferred parameters. 

Our results are in excellent agreement with X-ray mea- 
surements by Schmidt & Allen (2007), who found — 



320_j;^o 1^P<^ ^'^d C200 = 5.1I;^;g; this again suggests there 
are no significant systematic errors in our weak lensing 
analysis. The concentration C200 — 4.71^ 2 consistent 
with the CDM expectation of 4.1 ±0.4 ()Neto et al.ll2007D 
for a "relaxed" halo of this mass. 

3. HST IMAGING, SURFACE PHOTOMETRY, AND 
STRONG LENSING ANALYSIS 

3.1. Observations and Data Reduction 

AbeU 611 was observed by HST/ ACS with a 36 min 
integration in the F606W filter (ID 9270, PI: Allen). We 
obtained the data from the HST archive in reduced and 
drizzled form; however, the background level was not 
exactly uniform among chips. We thus remeasured and 
subtracted the background in the two relevant chips, tak- 
ing care to avoid bias from the extended halos of galaxies, 
particularly the BCG. 

3.2. BCG Surface Brightness Profile 

The distribution of stellar mass in the BCG must be 
measured to model strongly lensed features and stellar 
velocity dispersions. The BCG surface brightness pro- 
file from HST imaging was fit to an elliptic al R^^^ pro- 
file using the two-dimensional Galfit code (|Peng et al.l 
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I2002D . Cluster galaxies were excluded from the fit by 
enlarging their SExtractor ellipses by a factor of 5; 
extended arcs and diffraction spikes were masked manu- 
ally. The P SF was accounted for using a TinyTim model 
dKrisl^'iggS), and variance was estimated from standard 
CCD parameters scaled by the exposure map. The fit 
parameters are given in Table [3) The uncertainties are 
dominated by systematic effects, primarily the sky level 
and the fitting region. Both were varied to obtain the 
estimates in Tabled 

We use the i?^/^ fit in our dynamics modeling, and Fig- 
ure [5] shows that this indeed provides a good description 
of the data in the '^3-20 kpc region over which we mea- 
sure velocity dispersions. For strong lensing purposes, 
we adopt the dual pseudoisothermal elliptical mass dis- 
tribution (dPIE)^^, which is characterized by two length 
scales Vcore and rcut, with Vcore < fcut- The dPIE model 
in Table [3] was obtained by adopting the ellipticity pa- 
rameters from Galfit and fitting the surface brightness 
as a function of elliptical radius. As shown in Figure [5l 
this is a excellent fit over the entire strong lensing regime. 
We determine the total luminosity from the dPIE model. 
Galactic extinction was corrected using the E(B — V) va l- 
ues and extinction coefficients of iSchlegel et al.l (|1998( ). 
Tranforming to the Vega system, and using the trans- 
formation AmB- F606W = 0.65 ± 0.05 calculated by 
iTreu et al.l (fT999f) for a z = 0.292 passive galaxy, we find 
Me.vega = -23.89 ±0.15 and Lb = (5.6 ±0.8) x 10" Lq. 

3.3. Multiple Image Constraints 

The mass distribution from ~30-90 kpc is constrained 
by three multiply-imaged sources identified in Figure [H 
Spectroscopic redshifts are available for sources 1 and 2 
(Richard et al. 2009, in preparation). We conservatively 
adopted as constraints only the multiple images we con- 
sider most secure. However, objects are detected at all 
locations of counterimages predicted by the model, pro- 
viding a vital check. 

• Source 1. Five images are identified by a complex 
morphology with five distinguishable knots per im- 
age. We chose to include only one knot (knot A 
of Richard et al. 2009, in preparation) as a con- 
straint, both to avoid overweighting these images 
and for computational efficiency. We have verified 
that including all knots does not significantly af- 
fect our results. Images 1.4 and 1.5 are split by the 
galaxy they bracket, with a likely radial counterim- 
age indicated at the bottom of the inset in FigurelHl 
Although this image is predicted by our model, the 
predicted position is too uncertain to make a secure 
identification, and we therefore decline to include 
this image as a constraint. 

• Source 2. This source forms a giant tangential arc. 
Three merging images are identified from bright- 
ness peaks. 

• Source 3. Images 3.1 and 3.2 are identified by 
their common and distinct "horseshoe" morphol- 
ogy. The redshift has not been measured spectro- 
scopically, and the faintness of the images, along 

11 Also called truncated PIEMD (see Eb'asdottir et al. 2007) and 
distinct from the original model of .Kassiola fc Kovncr, I1993i) . 



with contamination from other sources, precludes 
accurate color measurements in our ground-based 
data. Three additional images are predicted by our 
initial models; two of these (3.3 and 3.4) have pre- 
dicted positions and fluxes that match well with 
features in the HST imaging, and we do include 
these as constraints. The third predicted image is 
in the cluster core, and although we identify a po- 
tential match in the inset of Figure [6l we do not 
consider the interpretation reliable enough for use 
as a constraint. 

In total we use ni — 12 images of =3 sources, 
resulting in 2{ni — Us) — 18 constraints on our model. 

3.4. Cluster Member Calaxies 

Cluster galaxies other than the BCG can perturb the 
positions of nearby images. These are included in our 
strong lens model via dPIE proffles, with properties 
set by scaling relations that are motivated by the fun- 
damental plane. The cluster member catalog was se- 
lected from resolved HST sources, detected using SEx- 
tractor, that have (1) projected cluster-centric ra- 
dius < 400 kpc (imaging is complete to ~200 kpc), (2) 
l^phot - Zciustorl < 0.2 (~2(t), (3) and mpeoew < 23. 
This selects 90 galaxies; an additional fainter perturber 
near image 1.3 was also included. SExtractor cen- 
troids and ellipticities were used, while the dPIE pa- 
rameters rcut and ctq were set by scaling relations. We 
ad opted the r*-band rela tion L/L^ — (cro/fo*)^'^^ found 
by iBernardi et al.l (|2003[ ) for early- type SDSS galaxies, 
where M^, = —21.39 and (Tq* = 158 km s~^ at the cluster 
redshift, accounting for luminosity evolution. However, 
recognizing the ~25% intrinsic scatter in this relation 
and possible systematic errors (e.g., photometry errors 
within the BCG), ag* was allowed to vary by imposing 
a Gaussian prior with a 50 km s~^ dispersion. We 
note that our limiting magnitude of mpeoew = 23 cor- 
responds to CT w 70 km s^^, and hence to a deflection 
angle of ~ 0'.'07 in the SIS approximation, well below 
our positional uncertainty of Upos = 0'.'5. 

The halo size of cluster members is uncertain, par- 
ticularly in cluster cores where tidal stripping may be 
very efficient. Natarajan et al. (2009) found that rc„t,* — 
45 ±5 kpc (3o-) in the core of CI 0024+16 from studies of 
galaxy-galaxy lensing. We therefore use the scaling re- 
lation L/L* = (rcut/rcut,*)'^, but impose a wide uniform 
prior on rcut,* (30-60 kpc). Similarly, rcore is taken to 
scale as L^^^, but rcore,* is fixed at 0.15 kpc, as this has 
a small effect. 

Several cluster galaxies are quite close to images used 
as strong lensing constraints. It is natural to wonder 
whether the fit may be improved by modeling them sep- 
arately. We have investigated freeing the 8 perturbers 
outlined in Figure [6] from the scaling laws by allowing 
their ctq and rcut parameters to vary separately (one per- 
turber at a time), and in no case were the predicted image 
positions significantly improved. 

We therefore judge the scaling relations sufficient to 
model the cluster galaxy population in the fit to strong 

1^ This assumes that the dPIE parameter cro is a good proxy for 
the measured velocity dispersion, as projected a nd averaged over 
an ap erture, which is true for our small rcore (see lEKasdottir et al.l 
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TABLE 3 

BCG Surface Brightness Profile 



Model 


Re or rcut 


rcore »TtF606W 


b/a 


PA 


ill/4 


9.79 ± 1.3 


16.77 ±0.12 


0.720 ± 0.008 


42.7 ±0.6 


dPIE 


10.0 ±0.8 


0.25 ±0.01 16.88 ±0.14 







Note. — Distances are in arcseconds. The scale is 4.33 kpc arcsec ^. The 
Re parameter from Galfit has been circularized by multiplying by (6/a)i/^. 
Uncertainties are primarily systematic and were estimated by varying the sky 
level and fitting region. 




Fig. 6. — HST/ ACS imaging of the cluster core (scaled logarithmically) and multiple image interpretation. Three multiply-imaged sources 
are identified by the first number in each image label, with decimals distinguishing separate images. Red circles have radii indicating the 
positional error from the best fitting model, crpoB = O'.'S. Cluster galaxies outli ned in yellow are potentially significant lensing perturbers, 
and the effect of freeing them from the scaling relations was investigated (see il3.4| l: this improved the fit only for the galaxy marked PI. 
Inset: BCG core with the R^^'^ model subtracted to reveal two additional likely radial counterimages used to verify the model. 



lensing data. However, when tension from kinematic 
data is introduced, the giant arc is better fit by free- 
ing perturber PI (Figure [6]) from the scahng laws. For 
consistency, we optimize PI in the strong lensing fit as 
well. For this perturb er, rcut is fixed by scaling from 
rcut,* = 45 kpc (Na taraian et al.ll2009r). and an is al- 
lowed to vary based on the iBernardi et al.l (|2003f l rela- 
tion, with a 50 km dispersion reflecting both the 
intrinsic scatter and the additional uncertainty arising 
from fixing rcut ■ It is not necessary to free both Vcut and 
do, since their effects are degenerate. 

3.5. Strong Lens Modeling of Abell 611 

Although image positions can in principle be deter- 
mined very precisely from HST imaging, modeling er- 
rors dominate in practice. These include unmodeled sub- 
structure and use of functional forms (e.g., (g)NFW) that 
may be inappropriate to describe the halo. The uncer- 
tainity CTpos in image positions must therefore be inflated 



from the astrometric error. This uncertainty also con- 
trols the relative weight of strong lensing in the combined 
analysis. We find our best-fitting models are unable to 
reproduce image positions to better than 0'.'37, and have 
chosen to adopt dpos = 0'.'5. Smaller values of a-pos would 
only strengthen our claim that NFW cannot simultane- 
ously fit all of our mass probes. 

Our mass model consists of several components (rele- 
vant parameters follow in parentheses): a pseudoellipti- 
cal (g)NFW model of the DM halo (e, 0, r^, cto.dm, /3), a 
dPIE model of the stellar mass in the BCG (cto.bcg), and 
dPIE models of the cluster galaxy perturbers described 
in §3.41 (r„,,f ^. an^. gn pi ). S ince the BCG structural pa- 
rameters were measured in §3.21 the only free parameter 
is o'o^BCG; or cquivalcnt ly Mi,/L = 1.57r g'Q gggrg^t 

/GL 

(iSand et a l. 2008; Elfas dottir et al.ll2007l ). Note that we 
model only the stellar mass in the BCG, and thus iden- 
tify the DM halo of the B CG with that of the cluster 
fsee lMiralda-Escu"dill995f) . Our models of perturbing 
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TABLE 4 

(g)NFW Fits to Strong Lensing Data 



Parameter 


Units 


Prior 


NFW Posterior 


gNFW Posterior 


(g)NFW DM halo 
e 

PA deg 

Ts kpc 
o-Q.DM km s"^ 
P 

M200 10" M0 
C200 


[0.1,0.3] 
[40, 48] 
[50, 800] 
[1000, 2200] 
[0.05,1.5] 


42.7 ±0.7 

136t?* 
1 256 + 30 

3.50+Si? 
10.0 ± 1.1 


0.220 ±0.018 
42.7 ±0.7 
68111 

n 44+0-40 

3 06+0-60 


dFlhl model 
f^O,BCG 


of BCG stars 
kms-i [109,345] 
solar 


< 216 

< 4.0 


< 239 

< 4.8 


Cluster galaxy perturbers 

o-o. kms-i (At, cr) = (158, 50) 
rent,, kpc [30,60] 
CTO.Pl kms-l (^t,cr) = (112,50) 


185 ± 18 
< 43 


179 ± 19 

< 45 
1461^1 


Source 3 redshift [1.5,2.5] 
Minimum XgL 

Number of constraints (10 parameters) 
Evidence ratio 


■^■^^-0.14, 

9.5 
1 


2.18 ± 0.17 
6.7 

18 

1.1 ±0.4 



Note. — Brackets indicate a uniform prior, while (^, a) indicates a Gaussian 
prior. Quantities without priors are derived, not directly inferred. Priors on e, PA, 
(To DM, and the source 3 redshift were determined from initial test runs; that on 
(To BCG corresponds to 1 < M,/Lb < 10. The inferred pseudoellipticity e implies 
b/a Ki 0.6 for the mass surface density. Limits on Mf/Lg reflect the dominant 
uncertainty in Mt, only, not in Lg. Error bars and upper limits refer to the 68% 
CL. 



cluster members, on the other hand, represent the total 
mass - dark and baryonic - since there is good evidence 
that th ese galaxies maintain individual halos as they fall 
inward (|Nataraian et al.|[200al . We assume the DM halo 
is centered on the BCG centroid, based on our weak lens- 
ing results and the coincidence of the Chandra (ID: 3194, 
PI: Allen) X-ray peak within the ~1ct astrometric uncer- 
tainty. Freeing the DM center in our strong-lensing fits 
results in a subarcsecond offset of < 3 kpc but reduced 
Bayesian evidence. The unknown redshift of source 3 
is also inferred with a uniform prior [1.5,2.5] based on 
initial test runs. 

We now have 11 free parameters, of which 4 are con- 
strained by physically motivated priors, and 18 con- 
straints. For each set of model parameters, Lenstool 
computes the image plane position r* of each image i, 
and then 



XSL = 51 



odd 



(5) 



Only image positions, not their fluxes or shapes, are used 
as constraints. 

The (g)NFW mass models inferred from strong lens- 
ing constraints only are given in Table |4l These one- 
dimensional marginalizations must be interpreted with 
caution, since severe correlations exist. These will be ex- 
plored in Sj5l where the results from our three mass probes 
are compared. We find that both NFW and gNFW mod- 
els fit the strong lensing data acceptably, with = 9.5 
and 6.7, respectively, for the maximum likelihood mod- 
els, corresponding to rms position errors of 0'.'37 and 
0'.'44. A shallow inner slope (3 — 0.44^q'24 is marginally 
preferred, but the additional complexity of the gNFW 



models is not demanded by the strong lensing data alone, 
as reflected in the evidence ratio of unity. 



3.6. Strong and Weak Lensing Combined 

At this point we can assess whether the constraints 
from strong and weak lensing are compatible. There are 
two reasons for such a comparison. First, we compare the 
mass enclosed within 100 kpc, approximately the bound- 
ary between the strong and weak lensing data, to deter- 
mine whether the lensing data are consistent where they 
(nearly) overlap. Since M3d,sl(< 100 kpc) = 3.041°;°^ x 

10^3 Mq and M3d,wl(< 100 kpc) 2.3t^J x lO^^ Mq 
differ by < Ict, we see no sign of a significant discrepancy. 

Second, we compare the mass enclosed within 1.5 Mpc, 
approximately the virial radius inferred from weak lens- 
ing, to determine if the (g)NFW profiles adequately de- 
scribe both lensing regimes simultaneously. This is di- 
rectly constrained from weak lensing to be M3d,wl(< 
1.5 Mpc) = (7.4 ± 1.4) X 10" Mq, and is extrapo- 
lated from the strong lensing results to be M3d.sl(< 
1.5 Mpc) = (3.70 ± 0.41) X lO^M© using an NFW pro- 
file. (Adopting a gNFW model instead shifts these by 
< la.) The significance of this difference 2.5cr) sug- 
gests marginal tension in the ability of a (g)NFW profile 
to fit all the lensing data simultaneously, but is not con- 
clusive in itself. 

The lensing data may be formally combined by setting 
= XsL + XwL- expected, little information on the 
inner slope is obtained: the posterior probability density 
is flat over a wide range {j3 ~ 0.6—1.5) and falls off slowly. 
Correspondingly, the evidence ratio NFW : gNFW = 
(2.2 ± 0.6) : 1 is near unity. 
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Radius [li,„ ' kpc] 




6400 6600 6800 7000 

Wavelength [A] (top: rest, bottom: observed) 

Fig. 7. — Central spatial bin of the BCG spectrum (black) com- 
pared with the best- fitting model described in the text (red) , with 
residuals shown below. The atmospheric B-band around 6900 A is 
excluded from the fit, as are narrow bright sky lines. Also masked 
are two poo r ly fit spectral lines, including Mg b. As discussed by 
IBarth et al.l l|2002l '). it is not possible to simultaneously fit Mg b 
and the Fe lines, since [Mg/Fe] is itself a strong function of the 
velocity dispersion. 

4. STELLAR KINEMATICS AND DYNAMICAL 
ANALYSIS 

4.1. Observations and Data Reduction 

We observed the Abell 611 BCG on 3 March 2008 us- 
ing the long-sHt mode on the LRIS spectrograph at the 
Keck I telescope, usmg a 600 mm^^ grating blazed at 
7500 A, targeting the redshifted Fe hnes around 6800 A. 
The 1'.'5 wide slit was ahgned within 3° of the major axis, 
as measured in i i3.2l The average seeing was ^1'.'4. Five 
exposures were taken for 130 min total integration. 

The LRIS data were reduced in a standard manner 
with bias subtraction, flat-fielding via quartz lamp expo- 
sures taken between science observations . and c osmic-ray 
rejection using LACOSMIC (jvan Dokkumi 120011 ). Spatial 
and spectral distortions were removed with the IRAF 

tasks IDENTIFY, REIDENTIFY, FITCOORDS, and TRANS- 
FORM. Wavelength calibration was determined directly 
from night sky lines, with residuals of order 0.2 A. The 
slit function was determined from night sky lines, which 
were then subtracted by low-order fitting. The spectral 
resolution was measured from unblended night sky lines 
to be cr = 99 km s"^ over the spectral range used for ve- 
locity dispersion measurement (6300-7100 A). The spec- 
tra were spatially binned to provide sufficient S /N (> 20 
per pixel) for a reliable determination. Bins were also re- 
quired to contain at least 5 CCD rows, comparable to the 
seeing element, to avoid an excess of correlated points. 
One bin was excluded due to possible contamination from 
an interloping galaxy. 

4.2. Velocity Dispersion Measurements 

Velocity dispersions were me asured using a direct fit - 
ting procedure similar to that of lSand et al.l ()2004l2008D . 
High-resolution template spectra of G0III-K2III giants, 
previously obtained using ESI on the Keck II telescope, 




Fig. 8. — Measured stellar velocity dispersions (rectangles) with 
NFW (lines) and gNFW (shaded) fits to strong lensing data alone 
(upper) and to strong lensing and kinematic data (lower). A NFW 
extrapolation from larger scales overpredicts the velocity disper- 
sion. Models have been PSF-convolved to match the data. For the 
data, radii have been circularized from their major axis positions 
by multiplying by (b/a)^^^. The dotted and dashed rectangles 
distinguish measurements on either side of the center. The 68% 
confidence regions are shown. 



were convolved with a Gaussian to match the LRIS in- 
strumental resolution. The templates were then red- 
shifted and rebinned to match the range (6300-7100 A) 
and dispersion (57 km s^^ pixel"^) of the BCG data. 
Model galaxy spectra were computed assuming a Gaus- 
sian line-of-sight velocity distribution (LOSVD). We fit 
these models to t he data with a modified version of the 
code described bv lvan der Marell(|1994 ). The best-fitting 
model is that which minimizes 



E 



S., - [Pa/ (A,) • (B,,„^,, o r)(AO + Pjv(A.)] 



(6) 

where Si is the measured spectrum, i enumerates the pix- 
els, Bai^^^z oT is the broadened and redshifted template, 
Pm{X) and P/v(A) are M and N order polynomials, re- 
spectively, and Ui is the flux uncertainty computed from 
standard CCD parameters. The free parameters are aios, 
redshift, and coefficients of Pjv and Pm- Pn is an addi- 
tive term used to model the galaxy continuum. Pm is a 
multiplicative term containing the normalization and any 
wavelength-dependent differences in instrument response 
between the galaxy and template spectra. 

Given the spectral range of 800 A, we anticipate w 9 
is necessary to model continuum variations on 100 A 
scales. Low N will cause a poor fit and a biased measure- 
ment, while high N allow the polynomial to fit spectral 
lines. Indeed, A^ = 9 provides a significantly better fit 
than N — 8, while for higher N the fit quality plateaus. 
We therefore take A^ = 9 and adopt a multiplicative or- 
der M = 2 to model any uncorrected differences in in- 
strument response. The best fit was obtained with the 
G9III template, and the central spatial bin is shown in 
Figure [71 Spectral regions around bright sky lines, the 
atmospheric B-band absorption, and two poorly-fit lines 
were excluded. The fit is excellent, particularly around 
the dominant Fe A5270 line, and the small, pattern-free 
residuals indicate no serious systematic effects. 
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TABLE 5 
BCG Stellar Velocity Dispersions 



Spatial bin 


(kpc) 


(S/N) /pixel 


crios (km s i) 


-5'.'2 to -3'.'7 


-22 to -16 


21 


316 ±67 


-3'.'7 to -2'.'6 


-16 to -11 


27 


308 ± 40 


-2'.'6 to -I'.'e 


-11 to -4.2 


39 


330 ± 28 


-I'.'e to -0('5 


-4.2 to -2.2 


58 


333 ± 18 


-0^5 to O'.'S 


-2.2 to 2.2 


73 


287 ± 15 


0^'5 to I'.'e 


2.2 to 4.2 


58 


293 ± 17 


to 2'.'6 


4.2 to 11.3 


38 


308 ± 25 


2^6 to 3'.' 7 


11 to 16 


25 




3^' 7 to 5'.' 4 


16 to 23 


21 


344 ± 40 


Note. — Nei 


gative (positi 


ive) positions 


are southwest 



(northeast) of the centroid along the slit. Errors are Icr 
statistical plus a systematic estimate (10 km s~^) added in 
quadrature. The 11-16 kpc bin is excluded due to possi- 
ble contamination from an interloper. The dispersion is 1.3 
A pixel" ^ . 

Measured dispersions for all bins are given in Table O 
The data are statistically consistent with a flat dispersion 
profile out to our radial limit. Rotation was not detected 
and must be < 20 km s~^, constituting < 0.5% of the 
kinetic energy. 

To evaluate systematic errors in our fitting procedure 
(due, e.g., to template mismatch or poor continuum fit- 
ting) , four different spectral regions were fit using various 
continuum orders and stellar templates (G7in-K0III), 
and the systematic uncertainty was estimated from the 
scatter in derived velocity dispersions to be ~ 10 km s~^. 
We add this in quadrature to the statistical errors derived 
from the surface to produce the estimates in Table El 
The effect of systematic changes to the velocity disper- 
sions on our results is discussed in t j5.2l 

4.3. Dynamical Modeling 

We have incorporated velocity dispersion constraints 
into Lenstool. Each model suggested by the MCMC 
sampler specifies the DM halo parameters and the stellar 
M^,/L. Projected stellar velocity dispersions are then 
computed based on the (g)NFW halo and the i?^/^ model 
of the BCG (i i3.2p . Spherical symmetry and isotropic or- 
bits are assumed here, as discussed extensively in i j5.2l 
and we therefore circularize the bin radii in our data. 
The dispersion and deprojected density profiles for the 
law are given bv lYound (fT976h . The effects of see- 
ing, slit widt h, and spat i al bin ning were accounted for as 
described in lSand et al.l (|2004i) . and 

2 X ^ / ^^i,obs model \ /rrx 
XVD = E { j 

was computed by summing over the 8 bins in Tabled 

Although these kinematic data alone cannot uniquely 
constrain the mass model, they acquire great power to 
determine the mass profile at < 20 kpc scales when 
combined with the strong lensing data. We coupled 
the kinematic and strong lensing data by setting ^ — 
XsL + XvD- total, this provides 26 constraints on the 
same (g)NFW models considered in tj3.5l The inferred 
parameters are given in Table[Sl Figure[5]presents the ve- 
locity dispersion profiles derived from NFW and gNFW 
models to these data, as well as those inferred from strong 



lensing data alone ( iJ3.5p . 

There are several important points to note. First, 
an NFW extrapolation of strong-lensing constraints on 
~30-90 kpc scales overestimates the mass at radii < 
20 kpc, and hence the central velocity dispersion. Sec- 
ond, an NFW fit can more nearly match the kinematic 
data when it too is imposed as a constraint, but the 
quality of fit at larger radii then suffers dramatically 
(AxsL — 22). The constraints on both scales cannot 
be simultaneously met. Third, when the asymptotic in- 
ner slope /3 is freed by adopting the gNFW profile, a 
similar fit to the kinematic data is obtained, while in- 
curring far less discrepancy with the strong lensing data 
(AxgL = 5). This is achieved by selecting /3 < 0.22 
(68%), a mu ch shallowe r densi ty profile than NFW. As 
discussed in iSand etldl ()2004[ ) , only upper limits on /3 
are available for systems without a radial arc. Although 
we identified potential radial arcs in §3.3( they have not 
been used to constrain our models, due to their uncertain 
identification. We therefore expect our inference of (i to 
be one-sided. 

5. JOINT ANALYSIS OF THE DARK MATTER 
DISTRIBUTION: 3 KPC - 3.25 MFC 

In this section, we now realize our goal of combining 
the three observational techniques: weak lensing from 
150 kpc to 3 Mpc, strong lensing from ~30-90 kpc, and 
stellar kinematics from ^-^3-20 kpc, in order to probe the 
dark matter distribution, independently of the baryonic 
contribution. The wide dynamic range of our data pro- 
vides a stringent test of theory. We seek to accomplish 
two goals: (i) to verify or otherwise the validity of the 
NFW profile for Abell 611, (ii) in the specific case of the 
gNFW profile, to determine the likely range of the vari- 
able inner slope In achieving these aims, we will also 
demonstrate the value of our approach. Specifically, we 
demonstrate that although an acceptable fit to the NFW 
profile might be achieved with a subset of our data over a 
more limited radial extent, such an agreement is illusory 
when the broader range of data is taken into account. 

We begin by considering Figure [9^, which shows the 
mass profile obtained from NFW fits to our three data 
sets. The weak lensing data is well-fit by an NFW profile, 
with a concentration approriate to its mass and in excel- 
lent agreement with X-ray data, as discussed in i )2.5l By 
looking at large scales only, we would therefore confirm 
CDM predictions. On intermediate scales, we found that 
positions of multiply-imaged sources could also be repro- 
duced to good precision by models including an NFW 
halo. However, this fit diverges from that obtained via 
weak lensing (and X-ray) data on Mpc scales, as indi- 
cated in Figure [IK. 

The most serious discrepancy comes on < 20 kpc 
scales, where we find NFW extrapolations from the 
strong lensing regime imply velocity dispersions signif- 
icantly higher than are measured. The dispersion data 
can be partially matched only by incurring drastically 
increased errors in the predicted positions of strongly- 
lensed sources. (Since strong lensing constraints are the 
most precise, even strong deviations are not visible on 
the scale of Figure [51 The severity of the problem was 
discussed in ^^3.51 ) 

We have also considered a simple generalization of the 
NFW profile, motivated by controversy in the theoretical 
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TABLE 6 
(g)NFW Fits to Strong Lensing 



Kinematic Data 



Parameter 


Units 


Prior 


NFW Posterior 


gNFW Posterior 


(g)NFW DM halo 










e 




[0.1,0.3] 


0.165 ± 0.008 


0.179 ± 0.009 


PA 


deg 


[40, 48] 


43.0 ±0.8 


42.8 ±0.7 




kpc 


[50,800] 






O"0,DM 


km 


[1000, 2200] 


1408 ± 30 


1663153 






[0.05,1.5] 




< 0.22 


jV/200 


IQl-* Mq 




fi 1+0.6 

'^•-'--0.5 


0+0.3 
■^■"-1.0 


C2OO 






7.05 ± 0.42 




dPlE model of BCG 
O"0.BCG 


stars 
km s"-"^ 
solar 


[109, 345] 


< 123 

< 1.3 


1691^2 
2.3 ±0.7 



Cluster galaxy perturbers 

CTQ* km s^-"^ 

rcut,* kpc 
(TO PI km s^"*^ 



{^l,a■) = (158,50) 
[30, 60] 
cr) = (112,50) 



167 ± 15 
169 ± 16 



157 ± 17 



162 



+13 



-17 
+(1.16 



Source 3 redshift 

„2 /,,2 



[1.5,2.5] 



Minimum Xsl/Xvd 
Number of constraints (10 parameters) 
Evidence ratio 



1.91 ±0.11 2.03_„ 11 

31.8/51.0 11.3/41.8 
18/8 

1 (1.9 ± 0.6) X 10^ 



Note. — See notes to Table |4] S ince rcut,- 
measure of its (fiat) posterior. 



essentially unconstrained, we omit a 



(g)NFW Fits to Strong 



TABLE 7 

f Weak Lensing and Kinematic Data 



Parameter 


Units 


Prior 


NFW Posterior 


gNFW Posterior 


(g)NFW DM halo 
e 

PA deg 

rs kpc 
o"0,DM km s-l 
13 

M200 1014 Mq 
C200 


[0.1,0.3] 
[40,48] 
[50, 800] 
[1000, 2200] 
[0.05,1.5] 


0.163 ± 0.008 
42.9 ±0.8 
236^24 

1415 ± 30 

6 2+"'^ 
6.95 ± 0.41 


n 170+0.009 

'-'•J^'J-O.OOT 
42 6+°-^ 
-,-,+16 
'1-6 

1682+11 
< 0.30 

q 0+0.4 

•^•''-0.6 


dFIE model 
o"0,BCG 
M*/Lb 


of BCG stars 
km s~i 
solar 


[109, 345] 


< 123 

< 1.3 


170+20 

l'^-30 
9 7+0.7 
^•'-0.8 


Cluster galaxy perturbers 
CTQ* km s~l 
rcut,* kpc 
o-Q.Pi km s^i 


= (158,50) 
[30,60] 
(Ai,(T) = (112,50) 


164 ± 15 
170 ± 16 


145 ± 26 
155 ± 16 



[1.5,2.5] 
(^, cr) = (0.81,0.05) 



Source 3 redshift 
Shear calibration m 

Minimum Xsl/Xvd/Xwl 

Number of constraints (ll parameters) 

Evidence ratio 



1.883 
0.85j 



1 qq^"-i'' 
J^-^^-O.ll 

r. 07 + 0.05 

1-°'-0.06 



-0.09 
-+0.05 
-0.07 

31.9/50.8/13081.0 16.0/38.2/13086.7 
18/8/13090 
1 (2.2 ± 1.0) X 104 



Note. — See notes to Table |6] 



literature regarding the inner form of the DJVI distribu- 
tion. By freeing the inner slope f3, gNFW models have 
the extra freedom to decrease the central mass as the 
kinematic data require, while incurring less damage to 
the fit quality on larger scales. We found only marginal 
evidence for a flat slope /3 < 0.7 (68%) in fits to strong 
lensing data alone. Including kinematic data in the inner 
20 kpc forced a flat slope with high confidence: (3 < 0.2 
(68%). Finally, we have formally combined our three 
datasets by setting 



2 2 I 2 I 2 

X = XWL + XSL + XVD 



(8) 



and sampling with Lenstool. Posterior distributions 
are listed in Table [T] The effect is to slightly lower the 
evidence in favor of the gNFW profile, since the NFW fit 
to strong lensing + kinematic data, when extrapolated 
to large scales, agrees somewhat better with the weak 
lensing data than does the gNFW fit. (This can be seen 
in Figure m) 

The Bayesian evidence overwhelmingly favors the 
gNFW model over the NFW model, by a factor of 2 x iC. 
This corresponds to a preference significant at more than 
99%. Correspondingly, a shallow inner slope is strongly 
preferred, with /3 < 0.56 at the 95% CL (/3 < 0.65, 99% 
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Fig. 9. — (a) Enclosed projected 2D mass inferred from NFW fits to our tliree data sets, witli curves braclseting the 68% confidence 
region. Note that a single NFW profile is unable to match the dark matter distribution at all radii; instead, the NFW-based extrapolations 
of data are discrepant on kpc and Mpc scales. The filled box gives the location of the giant tangential arc. For reference, crosses denote 
weak-lensing circular aperture masses obtained from the statistic (Clowe et al."1998) (68%). Note that these points are correlated. The 
filled diamond denotes the X-ray mass at the virial radius (Schmidt & Allen 2007) (68%). The stellar contribution of the BCG is shown for 
the M«/Lb of Table [6] Lensing fits are fully elliptical, and the radius plotted here is the elliptical radius, (b) Same, except using gNFW 
models. Note that the additional flexibility of the gNFW profile permits less mass on < 10 kpc scales. 



CL). The inner slope is degenerate with the scale ra- 
dius and the M^/L ratio of the baryons, as shown in 
Figure [ini but we cannot reconcile (3 with an NFW-like 
slope if the scale radius is to be compatible with the data, 
and the M*/L ratio consistent with stellar evoluti on the- 
ory. A lthough recent numerical work (e.g., Navarr o et al.l 
120081 ) has suggested that slopes shallower than NFW may 
be found just beyond current limit of numerical reso- 
lution, the difference may be insufficient to explain the 
data. FigurefTTldemonstrates this by showing that on the 
scales we probe, the Einasto pr ofile, another functio nal 
form fit to A^-body results (e.g.. iNavarro "eran[200l . IS 
much more similar to the NFW form than either is to 
the gNFW models we infer. 

Although the gNFW profile provides a simple alterna- 
tive that our data support decisively over NFW, it too 
cannot be considered a fully satisfactory description of 



the data. Tension with the strong lensing data remains 
( §3.5p . and the observed flatness of the dispersion data 
(Figure [8|) is not reproduced. The shape of this profile 
is likely affected by the poorly-understood ways in which 
baryons shape DM in cluster cores, which we discuss in 

m 

Finally, motivated by some theoretical suggestions that 
the NFW profile may in fa ct better describ e the to- 
tal matter distribution (e.g.. iGao et al.l 120041 ) . we have 
constructed the total density profile by adding the dark 
and stellar components of our combined gNFW models. 
This total density profile is indeed intriguingly closer to 
NFW. However, we note that when comparing fits to our 
two-component gNFW+BCG model with fits to a single- 
component NFW model, the two-component model is 
still favored by a factor of 545 ±301 in evidence. 

Since our analysis includes a number of simplifying as- 




Fig. 10. — Left: Degeneracy between DM inner slope (3 and scale radius Ts- The slope is compatible with NFW only for scale radii 
much larger than allowed by the data. Contours indicate the 68%, 95%, and 99% confidence regions. MCMC samples are plotted. Right: 
Degeneracy between /3 and M^/Lb- Large M^/Lb ratios can mimic a shallow slope f): they imply a more massive B CG, and hence less 
DM in the central region in order to maintain the same total mass. Fundamental plane constraints at this redshift from lTreu fc Koopmanj 
1I2OOI) are shown (68%). 




Fig. 11. — Left: Logarithmic slope of the density profile inferred from (g)NFW fits to weak and stron g lensing and s tellar kinematics. 
For comparison, Einasto profiles with a = 0.129 and 0.219 are plottedj bracketing the range noted by Navarro et al.' II2004I ) in cluster 
simulations; here r_2 is fixed to the median inferred in NFW fits. [Navarro ct al. (2004) used such a profile to fit the inner regions of 
clusters in their A'^-body simulations. All confidence regions are 68%. Right: Density profiles inferred from (g)NFW fits to the combined 
data set. The BCG profile is also shown, using the M, /L ratios inferred in the gNFW model, along with the X-ray gas profile provided by 
S. Allen f i|5.3l l. For reference, the dashed line gives a gNFW extrapolation of the hot gas outside the range of the X-ray data. 



sumptions that could in principle bias the inferred den- 
sity profile, particularly (3, we have performed a battery 
of tests to evaluate known systematic effects, which are 
discussed in turn below. We find that none of these sys- 
tematic effects can seriously alter our basic results, and 
thus believe that our findings are robust. FigurefT^shows 
the distribution of [3 we infer from our combined fit to 
all data sets and summarizes the effects of the system- 
atics discussed below. In all cases, we retain /3 < 1 at 
> 98% confidence. Furthermore, the dominant system- 
atic is likely to be the exclusion of baryons from the N- 
body simulations to which we compare data. Although 
the theory is insufficiently developed to make good nu- 
merical estimates, it is generally believed that including 
baryons will widen the discrepancy we find with simula- 
tions (see Sjni). 

5.1. Triaxiality and Projection Effects 



Dark matter halos are expected to be triaxial. In the 
thin- lens approximation applicable to galaxy clusters, 
lensing depends only on the projected surface density, 
which we model in a fully elliptical fashion. However, 
unmodeled elongation along the line of sight {los) does 
affect dynamical inferences, as these depend on the 3D 
mass structure. The possible impact of triaxiality on our 
results can be estimated by using priors from cosmolog- 
ical simulations to infer the likely 3D shape of the DM 
halo. 

The halo pseudeoellipticity inferred from strong lens- 
ing (Table |4|) co rresponds to a pro.iec ted (mass) axis ra- 
tio oi q = 0.6 (jGolse fc KneibI 120021 Eqn. 28). Equa- 
tions relating the shape of a triaxial ellipsoid, with axis 
ratios c < & < 1, to its 2D p rojection are given by 
iRomanowskv fc Kochanekl ()1998D . We used Monte Carlo 
sampling to find the 6, c and los orientations compati- 
ble with the observed q and prior distributions p{c) and 
p{b\c) from the A^-body DM simulations presented by 
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DM inner slope P 

Fig. 12. — Marginalized posterior probability density for 0, as 
inferred from the combined fit to weak lensing, strong lensing, and 
kinemat ic data. The effects of the systematic effects discussed in 
^5.H5.4l on the 68% upper limit are illustrated by arrows. 

IJing fc Sutd ((200l . as SUES ested bv lGavazzil (|2OO50 . We 
find that the radius along the los is > 0.6 (2cr) of the aver- 
age radius in the plane of the sky; that is, strong compres- 
sion along the los is unlikely Furthermore, interaction 
with baryons tends to make halos roun der and less triax- 
ial th an those i n pure DM simulations ([Gustafsson et al.l 
[200l [D ebattis ta et al.ll2008l : lAbadi et al.ll2009f) . 

iGavazzit (|200l )" has estimated the error on the enclosed 
mass as inferred from kinematic data by assuming spher- 
ical symmetry, when the actual mass distribution is a 
prolate (oblate) ellipsoid aligned with major (minor) axis 
along the los. A prolate halo would exacerbate the dis- 
crepancy we find with NFW, since stars move faster 
along the major axis, and hence dynamical masses are 
overestimated. An oblate halo with short axis along the 
los causes an downward bias in dynamical mass, which 
is < 30% for the above limit on c. This bias can be ap- 
proximately compensated by increasing the velocity dis- 
persion measurements by 15%. Repeating our combined 
gNFW analysis with this modification yields (3 < 0.57 
(68%; (3 < 0.92, 95% CL). 

We therefore believe our results are robust against 
likely projection effects, although the significances may 
vary. Nevertheless, recognizing that triaxiality has a 
complicated effect on dynamical masses, and the pos- 
sible inadequacy of simple models in the present situa- 
tion (e.g., the simplifying assumptions bv iGavazzil 120051 
of a spheroidal halo aligned with the los and massless 
tracers), we intend to pursue this issue more rigorously 
in future work, both observationally and in modeling. 
Additionally, a larger sample will allow us to assess the 
likelihood that chance alignments of triaxial halos can 
explain our findings. 

5.2. Velocity Dispersion Measurements and Modeling 

In our dynamical modeling, we assumed isotropic stel- 
lar orbits and a Gaussian LOSVD. Since the mass distri- 
bution inferred from a velocity dispersion profile depends 
on the anisotropy tensor, the effects of these assumptions 
must be considered. Orbital structure has been stud- 
ied extensively in local samples of cD, clus t er, and field 
ellipt ic al g alaxies f e.g. iGerhard et al l 1998; 'Saglia et al.l 
120001: iKronawitter et all 120001 : IGerhard ct al.. .2001, : 



iThomas et alll2005l[2007l see also lKelson et al.ll2002l and 

references), and also in a few distant early- type galaxies 
(|Treu fc Koopmansl[200l . The broad consensus is that 
orbital structure in the inner regions is remarkably con- 
sistent: along the major axis, orbits range from isotropic 
to slightly radially biased, with the isotropy parameter 
f3 typically < 0.3 and up to « 0.5. Moreover, devia- 
tions from a Gaussian LOSVD are small, with /13 and 
typically a few percent and (rarely) up to 10%. At 
the radii we probe (< 0.4i?e), radial bias causes the pro- 
jected velocity dispersion to be higher than would be 
seen for isotropic orbits in the same mass distribution 
(|Binnev fc Tremainel 119871 . Fig. 4-13). Accounting for 
radial anisotropy would thus strengthen our results (i.e., 
force (3 down ward). Tangential b ias at these radii is quite 
rare; however. ISand et al.l (|2004f ) considered this possibil- 
ity and estimated (3 would be increased by < 0.2 in the 
clusters they studied. We conclude that our results are 
robust to any plausible level of orbital anisotropy. 

In i)4.21 possible systematic errors stemming from tem- 
plate mismatch and uncertain continuum fitting were 
found for our dispersion measurements at the 3% level. 
Additionally, we note a possible systematic difference of 
up to '^10% for measurements on either side of the BCG 
center. Since our basic findings were unchanged when 
the dispersions were increased by 15% ( ij5.ip . we do not 
believe our results are substantially impacted by these 
uncertainties. 

5.3. Intracluster Gas 

We have neglected the intracluster gas in our mass 
models, and our measured dark matter pr ofiles there- 
fore include a small contribution from it. lAUen et al.l 
(|2068l) concluded from X-ray studies that the average 
gas mass fraction within r25oo ~ 500 kpc was 0.10 ± 0.01 
in Abell 611, which decreases toward smaller radii. To 
assess the impact of neglecting gas on < 20 kpc scales, we 
subtracted the gas profile (kindly provided by S. Allen) 
from our gNFW models and found the slope was un- 
changed (A/3 ~ 0.01), as expected based on the very 
low gas fraction at small radii (Figure [TT]) . We caution 
that this is based on an inward extrapolation of the gas 
profile, which is not measured for r < 17 kpc. 

5.4. Other Observational and Modeling Errors 

We have verified that our results on /3 are not overly 
sensitive to the measured BCG size (i?e and rcut) by per- 
turbing it by the systematic error estimated in Table [31 
We have also made gross changes to seeing by perturbing 
it by 30%. In all cases the upper 68% confidence limit 
on /3 increased < 0.04. Finally, we note that although 
introducing ellipticity via the potential can produce un- 
physical negative surface densities at large radii (i.e., the 
outer regions o f our wea k lensing da ta) when Ts is suf- 
ficiently smaU (jGolse fc Kneib 2002.: ISand et al.l [2008.) , 
we checked that the reduced shear nevertheless remains 
physical in our models. 

6. DISCUSSION 

Our analysis of the mass density profile of Abell 611 
over an unprecedented range in radii gives two funda- 
mental results: an NFW profile is inconsistent with the 
data, and the logarithmic slope of the inner dark matter 



16 



Newman et al. 



density profile is < 0.3 (68%). It should be noted that our 
methodology allows us to disentangle the relative contri- 
butions of dark matter and baryons to the mass density 
profile, and our measurements therefore refer truly to 
the dark matter component. However, the cosmological 
A^-body simulations to which we compare our measure- 
ments do not include baryons. Although their impact is 
thought to be smaller in galaxy clusters than in galax- 
ies, baryons can be gravitationally dominant in cluster 
cores and may thus alter the dark matter distribution, 
particularly the inner slope. 

High-resolution, cosmological A^- body -I- gas dynam- 
ical simulatio ns of galaxy clusters (iGnedin et al.l l2004f ) 



and gal axies (lGustafsson et al.1 120061 : lAbadi et aL 1 120091 . 
but see iRomano-Diaz et al.l 20081) have shown that as 
baryons condense in the center, they deepen the poten- 
tial and steepen the dark matter distribution. In this 
scenario, the discrepancy we find with NFW-like dark 
matter profiles would be enhanced. 

IGnedin et al.l (|2004f ) confirmed that this steepening oc- 
curs whether or not baryon cooling (as well as star forma- 
tion, supernova feedback, and heating by the extragalac- 
tic UV background) was included, although the effect 
is larger when cooling is allowed. Cooling is certainly 
very important in cluster cores, but modeling the degree 
of cooling is currently limited by poor knowledge of a 
variety of baryonic processes, including viscosity (e.g., 
IDolag et a l. 2005; Puchwcin ct al. 2005) and AGN feed- 
back fe.g.. IPuchwein et al.l 120081) . IGnedin et all (|2004D 
acknowledge the common "overcooling problem" in their 
simulations, which should therefore provide an upper 
limit to the steepening. They found that a simple 
model of adia batic contraction, modifie d from the orig- 
inal model of IBlumenthal et all (1986") , reasonably de- 
scribes the behavior of dark matter in their simulations. 
This has been implemented in the public Contra code, 
which we used to estimate the response of the dark mat- 
ter halo in Abell 611 to the BCG. We find the dark mat- 
ter mass enclosed within 10 kpc, approximately the in- 
ner limit of the simulations, is enhanced by ~45%. Since 
the theoretical and modeling uncertainties are large, this 
clearly should be considered a rough figure. However, 
it does demonstrate that substantial contraction (steep- 
ening) of the dark matter is possible, and we note that 
this uncertainty may be comparable or larger than possi- 
ble oppositely-directed biases from unmodeled triaxiality 

(gED- 

Several researchers have suggested that other physi- 
cal processes, particularly dynamical friction between in- 
falling baryons and the host halo, may counter and sur- 
pass the effects of adiabatic contraction in galaxy clus- 
ters. This work has been either (semi)-analytic, or based 
on A/'-body simulations without hydrodynamics, cooling, 
or (typically) a cosmological context, but is neverthe- 
less important to bear in mind. If dynamical friction 
acts as strongly as claimed, our results could then be ex- 
plained with no contradiction to cosmological A'^-body 
simulations of dark matter: initial NFW halos would 
form flatter co res through their interaction with baryons. 
lEl-Zant et all 12001, 2004) considered an initially clumpy 
distribution of baryons and found they were efficient at 
"heating" the halo via dynamical friction, which softened 
the cusp. Their treatment of pure baryon clumps as un- 
strippable point masses may, however, exaggerate this 



effect. Indeed, iMa fc Bovlan-KolchinI (|2004f) found that 
infalling DM subhaloes, depending on their mass and 
con centration, could st eepen or soften the net DM cusp, 
and lNipoti et al.l (j2004D found a similar result depending 
on the baryon fraction of infalling galaxies. The mecha- 
nism by which the baryon clumps maintain their energy 
over sufficient timescales, without fragmenting and form- 
ing stars, is also unclear (jMashchenko et al]|2006 ) . As a 
way of understanding its "universalitv." lGao et all()2004[ ) 
suggested that the NFW profile is a dynamical attractor, 
to which coUisionless matter (stars and DM) is driven in 
a hierarchical formation picture. Since stars are centrally 
domi nant, this would imp ly a shallow DM profile. How- 
ever, [Gnedi^^Lari (|2004f) did not confirm this hypothesis 
in simulations including baryons. A recent analysis of a 
broad range of cosmological simulations (Gnedin 2009, 
private communication) - including a variety of condi- 
tions and detailed treatments of the relevant physics - 
showed that baryonic physics can increase the range of 
contraction, but does not change the sign, i.e. does not 
lead to "expansion." 

We thus conclude that our observations are hard to 
reconcile with the current understanding of cluster for- 
mation in the context of a ACDM cosmology. Either 
some of the assumptions about the dark matter back- 
bone are incorrect or - more likely in our opinion, given 
the number of independent lines of evidence supporting 
CDM on large scales - baryonic physics is not yet suf- 
ficiently well understood. Given the centrality of this 
problem and the success of CDM at other scales, this 
seems a worthy goal to pursue, although its solution will 
likely require advances in computing power, algorithms, 
and understanding of the relevant physics. 

A final caveat is that cluster mass density profiles may 
not be universal. Rather, a broader distribution of inner 
slopes may exist in nature than that predicted by pure 
dark matter simulations, potential ly depending on the 
merger history of the cluster (see iNavarro "etaI1l2008L 
and suggestions by our previous analysis of three clusters 
in Sand ct al. 2004). If that is the case, it is imperative 
to collect data for a larger sample of clusters, of similar 
quality to those presented here for Abell 611, in order 
to characterize the moments of the distribution as well 
as the mean. Ideally the next generation of simulations 
including baryons should be powerful enough to enable 
the simulation of large numbers of clusters and therefore 
allow a comparison with the measured distribution of 
inner slopes, taking into account selection effects. 

7. SUMMARY 

We have constrained the DM profile of Abell 611 from 
3 kpc to 3 Mpc and find that an NFW profile cannot 
simultaneously fit our lensing and kinematic data. We 
confirm the necessity of using mass probes over a wide 
range in cluster-centric radius to make the strongest com- 
parisons with theory. Freeing the DM inner slope f3 in 
our models increases the evidence by a factor of 2 x 10'* 
(i.e., at more than 99% confidence) and selects a shallow 
slope (3 < 0.3 (68%). We intend to apply our technique 
to a wider sample of clusters in order to test the univer- 
sality of our findings and measure the intrinsic scatter in 
cluster mass distributions. 
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Fig. 13. — (a) Stellar cUipticities sliowing large-scale, smooth variation of the PSF across the focal plane. Light regions correspond to the 
excluded CCD chip (upper left) and masked stellar halos. (b) Residual ellipticities, small and pattern-free, after subtracting the polynomial 
model described in the text. 



APPENDIX 
POINT SPREAD FUNCTION 

In order to measure the small shape distortions to background galaxies in the weak lensing regime, we must correct 
the effects of the point spread function (PSF), which is anisotropic and varies across the focal plane. We select bright 
but unsaturated stars from the stellar branch of the Rc^rh plane, where r/j is the half-light radius, and us e pietshapes 
from t he Imcat suite to measure ellipticities and polarization tensors. Our method is based on that of iKaiser et al.l 
(|1995l KSB). In this scheme, shape measurements are performed with a window function to avoid formally infinite 
noise properties. We choose a Gaussian window with cr = = r;i/(2 In 2)^/^, which is optimal in the case of a Gaussian 
source. For stellar measurements, Vg is fixed at the stellar median of (r*) = 1.8 pixels (0'.'36). For galaxies, this is the 
minimum allowed rg. For every object, the local (constant) background as estimated by SExtractor is subtracted. 
Ellipticities are then defined by Cq — {Qii — Qti^ Qi2}/('5ii + Q22), where Qap are the weighted quadrupole moments. 

Stars used for PSF measurement must have (1) tr P*^ > and tr P*^ > 0, where P*,„ and P*^ are the smear and 
shear polarizability tensors, respectively, (2) no masked pixels or pixels assigned by SExtractor to neighbors within 
2rg, (3) |e| < 0.2, and (4) d < 0.1 pixel, where d is the distance between the centroid computed with and without a 
weight function, as discussed in §2.31 

A two-dimensional polynomial of degree 7 is fit to e|, 62, and tr P*„ with iterative cr-clip rejection and used to 
interpolate the 422 stars throughout the field. The measured stellar ellipticity field and the polynomial fit residuals are 
shown in Figure [T5] Note that the large-scale, smooth PSF variations are removed. The PSF variance is reduced from 
iTea ~ 0.01 to « 5 X lO^'^. At the position of each object, the interpolated stellar anisotropy kernel g* is computed 
using the trace approximation to P,*™: 

= </(i tr P;„). (Al) 

The anisotropy-corrected ellipticity e' is then 
where Psm is assumed diagonal. 

GALAXY SHAPE MEASUREMENT 



In addition to the anisotropy introduced by the PSF and corrected with Equation I A21 galaxy ellipticities must be 
corrected for the isotropic sme aring caused by see i ng an d the window function. This is done using the preseeing shear 
polarizability P^ described bv lLuppino &: Kaisej ()1997l ). We use the trace approximation 

P''^l (^tr P,,.,„^ - tr P„„,„^ ^ trP*"'"^ )) ' ^^^^ 

where the quantity in angled brackets is the median stellar value. Since P'^ is a noisy estimator, we correct galaxy 
shapes using an average value over galaxies with similar properties. Specifically, P"' is fit to a three-dimensional 
polynomial that is quadratic in Pp, rh, and |e|, with iterative rejection. We find this is adequate to model the 
variation in P'^ . The interpolated P'^ is used to estimate the reduced shear = 7a/(l — k) by 

9c = e'jP\ (B2) 
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The uncertainty in is also estimated from like galaxies. For each galaxy, the nearest 50 neighbors in the Rc^rh 
plane are selected, where the distance along each axis is normalized by the standard deviation. The dispersion in ga 
among these neighbors is taken as the uncertainty. We note that this uncertainty (median cxg = 0.24) is dominated by 
the randomly distributed intrinsic shapes of galaxies. 
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